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THE  REFRACTION  OF  A PLANE  SHOCK  WAVE  AT  AN 
AIR-WATER  INTERFACE 
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INTRODUCTION 

The  nonnel  refraction  at  an  air-water  Interface  of  a plane  ahock 
wave  Incident  through  the  air  la  Inveatlgated.  The  purpoaea  of  thla 
Inveatlgatlon  are  twofold: 

a.  To  define  under  what  condltlona  the  In-water  ahock  front 
and  particle  velocltlea  and  the  changea  In  denalty  can 
be  predicted  by  acouatlc  theory;  and, 

b.  To  determine  under  what  condltlona  the  refracted  preaeure 
can  be  calculated  by  treating  the  Interface  aa  a rigid 
wall. 


Jump  Condltlona  Acroaa  A Shock 

the  Jump  condltlona  acroaa  a ahock  expraaa  the  conaervatlona  of 
maaa,  momentum,  and  energy  In  a control  volume  encloalng  the  ahock. 
From  a reference  frame  moving  with  the  ahock,  theae  condltlona  are 
given  by  [1] 

Oa^a  “ Ob^b  (la) 

2 

‘‘a  PaVa^  " ’’ b Ob  ^b  (lb) 

®a  /2  ■ eb 

Here,  the  aubacrlpta  "a"  and  "b"  refer  to  ahead  of  and  behind  the 
ahock,  reapectlvely.  The  aymbola  p,  P,  e,  and  v denote.  In  order, 
the  fluid  denalty,  preaeure,  apeclflc  Internal  energy,  and  relative 
velocity  with  reapect  to  the  ahock. 

The  former  three  of  theae  quantltlea  cannot  vary  Independently 
acroaa  the  ahock  but  are  constrained  to  aatlafy  the  equation  of  atate 
of  the  fluid.  For  air,  the  perfect  gaa  equation  of  atate  ylelda  the 
relatlonahlp 


e - P/((-r-l)  P I (2) 

where  > la  the  ratio  of  apeclflc  heata.  For  water,  the  conatralnt 
la  expreaaed  by  the  Talt  equation  of  atate  [2] 

P - Pq  + B ((  <*/  Pwo)"  - 11  (3a) 

Hare,  ^wo  la  the  water  denalty  at  atmoapherle  praaaura  Pq,  n la  a 
NoU:  Manuicript  MibmitUd  SapUmbar  80,  1977. 


1 


dlMnsionlcss  constant,  and  B la  a aacond  conatanC  daflnad  by 


B - ‘’wo 


(3b) 


vharc  ia  tha  apaad  of  sound  In  water  undar  atoMspharlc  conditions. 

Mechanical  Shock  Conditions 

Equation  (la)  glvea  the  condition 

Va  - 1)  Vb 

where  the  shock  compressibility  factor  1)  la  defined  by 


(4) 


^ • Pb^o, 

Upon  substitution  of  equation  (4),  equation  (lb)  ylelda 

V*b  " <*’a'0a>  (1/11(11-1)1 

where  the  preasure  Jump  factor  y la  given  by 

y - Pb^Pa 


(5) 


(6) 


(7) 


Equatlona  (4)  and  (6)  are  known  aa  the  mechanical  shock  conditions 
since  they  are  Independent  of  the  fluid  through  «rhlch  the  shock  la 
propagating.  Other  useful  relations  obtained  directly  from  these  con- 
ditions are 

2 


(Po/pa>(y-l)  [V(11-1)1 


(8a) 


(V.-Vb)=-  (P,/pa)  (y-1)  ((11-1)/T11 


(8b) 


To  proceed  further,  the  thermodynamic  propertlaa  of  the  medium 
through  which  the  shock  Is  moving  must  be  considered. 

Solution  For  Air 

Equation  (Ic),  which  expresses  conservation  of  energy  acroas  the 
shock,  becosMS,  on  utilising  the  thermodynamic  conatralnt  defined  by 
equation  (2), 

Iv/(vl)HPa/ps)  + v\/2-(v/(v-l)J(Pb/pb)  +Vb/2 


(9) 


j 


I 

\ 


I 

I 


I 

»■ 

it 

y 

i 
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On  substituting  equations  (6)  and  (8a)  and  dividing  through  by 
equation  (9)  yields 


. Ticy-n 

Y-1  ■^2(71-1) 


_y— 
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♦ y-l 
2T1(T1-1) 


(10) 


Some  simple  algebraic  manipulation  of  this  equation  gives  the  relation 
between  the  compressibility  factor  71  and  tha  pressure  Jump  factor 
y for  an  air  shock  as 


T|  “ (M-y)/(l  My) 


(11a) 


where 

M - (y-1)/(y  + 1) 


(11b) 


Solution  for  Water 

For  a water  shock,  the  relation  between  1)  and  y la  obtained 
directly  from  the  Talt  equation  of  state,  equation  (3a),  as 

^/n 

“H  - [(y  -»■  8)/(l  + 3)J  (I2a) 

where 

8 - (B  - Pq)  / P,  (12b) 

Shock  Refraction  at  an  Air-Water  Interface 


The  refraction  at  an  air-water  Interface  of  a plane  shock  wave 
Incident  through  the  air  la  schematically  Illustrated  In  Figure  1. 

From  Figure  1,  the  relative  velocities  of  the  fluids  as  seen  from 
the  shocks  can  readily  be  ascertained.  For  the  refracted  water  shock, 
these  relative  velocities  are: 


V — U" 

a 


(13a) 
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Iiirf«r  I 


The  water  shock  compressibility  factor  'Hw  is  found  from  equation  (I2a) 

as 

’Tw  - Dw2/  Owo-  [ 1 + n a(y  - 1)  j 


where 


or  - Pq/uB  - Po/PwoCwo 

• 

(18) 

and  where  equation  (12b)  has  been  used  to  obtain 

B - B/Po  - 1 - 1/n  cr  - 1 

(19) 

* 

Similarly,  for  the  incoming  air  shock,  the  shock  front  and  particle 
velocities  are  calculated  from  equations  (8)  and  (11a)  as 

U - j (Po/Pao)  [ (x-tM)/(l-M)j  j (20a) 

ui  - (x-1)  j(Po/pao)  [(l-M)/(x  +M)jj 

(20b) 

t 

< 

where 

t 

i 

• »■ 

X - Pi/Po 

(21) 

1 

For  the  refracted  air  shock,  the  relative  velocities  with 
to  the  shock  are  seen,  from  Figure  1,  to  be 

respect 

- 

V J ■ u J - U* 

(22a) 

j 

v^  -U2  " U' 

(22b) 

1 

1 

Hence, 

1 

v -V  ■ u,  _ u^ 
a b 1 2 

(23) 

j 

i 

i 

where  u^  and  u,  are  given  by  equations  (20b)  and  (15b),  respectively. 
However,  from  equations  (8b)  and  (11a),  ~ also  be  written  as 

• 

^a  - "'b  “ <*-!>  (Pl'Pai)  [(l-M)/(»-m)l 

r 

1 

(24a) 

i 

where 

* - ?2/?i  - (P2/Po)/(Pl/Po>"y/* 

(24b) 

Equation  (23)  than  becomes,  on  substituting  the  appropriate  expressions. 
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where  % Is  given  by  equation  (17).  Upon  multiplying  through  by 
( pao/Po  substituting  z"y/x,  end  noting  from  equation  (11a) 

that 


P.l/Pao  “ (M  + X)  / (1  + MX) 


(26) 


equation  (25)  yields 


(y-x)  fd-MXl-Wx)  1 -(x-W.(ljjD_\^  - 

[(y+Mt)(M4x)  J V*  + X / \Pwo  / 


(y-1)^ 


The  solution  of  equation  (27)  gives  the  refracted  pressure  ratio 
y (~P2/Po)  terms  of  the  incoming  pressure  ratio  x(*Pi/Po). 


Comparison  Between  The  In-Water  Shock  Waves  Obtained  Prom  Exact 
Shock  And  Acoustic  Theories 


The  process  of  shock  refraction  produces  at  the  water  surface  a 
step  change  in  pressure  with  a magnitude  given  by  P2  ' Pq  **  Pq  • 

The  response  of  the  water,  treated  as  a linear  acoustic  medium,  to 
this  step  change  in  surface  pressure  is  well  known  [3].  Retaining  the 
previous  notation  and  denoting  values  predicted  by  acoustic  theory 
with  a superscript  the  propagation  velocity  U"*  of  the  step  (shock 

front)  into  the  water  is  given  by 

U"*  - Cwo  (28) 

Behind  this  front,  the  water  particle  velocity  uf  and  the  water  den- 
sity p*w2  found,  respectively,  as 

u*2  * (Po/pwo  C|^)  (yl)  (29) 

and 

P^/  Pwo  • 1 + cr(y-l)  (30) 

with  a defined  by  equation  (18). 

The  exact  shock  and  linear  acoustic  predictions  of  the  water  res- 
ponse to  the  refracted  pressure  wave  ^ can  be  compared  by  forming 
the  ratios  U"/U*'*,  U2/U2t  snd  Pw2/p'*' *2  From  equations  (15), 

(17),  (28),  (29),  and  (30),  these  ratios  are  determined  as 


IT'/U"* 

. j Q (1  t 

1 <1  + nQ)l/n  - 1 J 

(3l«) 

U2/U*2 

. ( LI  ^ 

j Q (1  ilQ)‘'“  j 

(31b) 

”■'> 

where 

Q ■■  or  (y-l>  Old) 

The  above  ratios  are  plotted  In  Figure  2 as  a function  of  the  re- 
fracted pressure  ratio  y ■ P2^1*o*  obtaining  these  curves,  values 
of^n-T.lS  (ref.  f21)  an^  a - 4.5  x 10*5  . 1013  x 10^  dynes/ 

cm^,  .-^wo  ■ l.O  gro/cm' , Cwo  ■ 1.5  x 105  cm/sec)  were  used.  For 

refracted  pressures  up  to  1200  atmospheres,  Figure  2 shows  that  the 
In-water  shock  front  and  particle  velocities  calculated  by  exact  shock 
theory  differ  by  less  than  107.  from  the  In-water  shock  front  and 
particle  velocities  predicted  by  acoustic  theory.  Figure  2 also  shows 
that  the  changes  In  density  predicted  by  exact  shock  and  acoustic 
theories  differ  by  less  than  IX  for  refracted  pressures  up  to  1200 
atmospheres. 

Comparison  Between  the  Refracted  Pressure  and  the  Pressure  Reflected 
From  a Rigid  Wall 

For  an  Incoming  shock  of  strength  x^Pi/P^  (see  Figure  1>,  the 
pressure  P211  reflected  from  a rigid  wall  Is  given  by  [1]. 

yp  - •’2R^*’o  " ^ 1(1+2M)  x-M)/  (l4Mx)  (32) 

For  the  same  Incoming  shock  strength,  the  pressure  P2  refracted  at  the 
air-water  Interface  Is  found  as  the  solution  of  equation  (27). 

To  compare  the  refracted  pressure  to  the  pressure  reflected  at  a 
rigid  wall,  the  ratio  P2/P2R  formed.  This  ratio  is  plotted  In 
Figure  3 as  a function  of  the  reflected  shock  strength  P2r/Po'  (3^* 
numbers  In  parenthese  below  the  abscissa  correspond  to  the  incoming 
shock  strength  P|^/P  required  to  yield  the  reflected  shock  strength 
P2r/Po.)  In  obtalnfne  this  curve,  values  of  M"l/6  (v"l.4)  and 

■ 1.2  x 10*3  gm/cm^  were  used.  Figure  3 clearly  shows  that 
the  ratio  of  the  pressure  refracted  at  an  air-water  interface  to  the 
pressure  reflected  from  a rigid  wall  remains  above  967.  for  reflected 
pressure  up  to  1200  atmospheres  (Incoming  pressure  up  to  156  atmos- 
pheres) . 
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CONCLUSIONS 


In-water  shock  front  velocity  U,  particle  velocity  u and  density 
p can  be  predicted  by  acoustic  theory  up  to  1200  atm  with  a maximum 
error  of  lOZ  at  the  high  pressures. 

I 

The  refracted  pressure  can  be  calculated  by  treating  the  surface 
as  a rigid  wall  up  to  1200  atm  with  maximum  error  of  4Z  at  1200  atm. 
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Fig.  1 - Schematic  illustration  of  the  refraction  of  a plane  shock 
at  an  air-water  interface.  The  symbols  U and  u denote,  respective- 
ly, the  shock  front  and  particle  velocities  with  respect  to  a fixed 
frame. 
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Fig.  2 — Compamon  of  exact  ihock  (thia  paper)  to  acouatic  theory 
for  water  ahockwave  velocity  U ",  particle  velocity  U2  and  water  den- 
aity  * denotM  acouatic  valuea.  Comparative  valuea  are 

plotted aa  a function  of  the  retracted  preaaure  ration  P2/Po‘ 
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Fig.  3 — R«tio  of  lefiacted  pRi6uic  Pj 
to  the  pieesiue  reflected  at  a rigid  waU 
aa  a function  of  the  reflected  ahock 
•t**"****  PaR/**©- 


